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Use of HERA jet production cross-sections can extend our knowledge of the gluon 
parton distribution function and yield accurate measurements of as{Mz) in addi- 
tion to illustrating the running of as with scale. 



1. Introduction 

The QCD processes that give rise to scahng violations in DIS at leading 
order, namely QCD Compton (QCDC) and boson-gluon fusion (BGF), 
also give rise to distinct jets in the final state provided that the energy and 
momentum transfer are large enough. These processes are illustrated in 
Fig. 1. Both BGF and QCDC depend on as{Mz), but the BGF process also 
depends directly on the gluon parton distribution function (PDF), whereas 
QCDC depends on quark PDFs. The QCDC process dominates at large 
Q^, where the most important quark PDFs are the well-determined u and 
d valence distributions (for HERA kinematics). Thus measuring jet rates 





Figure 1. Leading order QCD diagrams from dijet production in DIS. (a) QCD Comp- 
ton; (b) Boson-gluon fusion. 
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at large allows a determination of a., (M^), with reduced uncertainty 
from the less well known gluon PDF. Such measurements of as{Mz) are 
described in Sec. 2. 

Conversely, if as{Mz) is well known, then the BGF process gives direct 
information on the gluon PDF. Historically, PDF fits have been made using 
NLO QCD in the DGLAP formalism to predict the inclusive cross-sections 
for world DIS data, but the gluon PDF affects these cross sections only 
indirectly through the scaling violations. Inputting jet data to such fits can 
improve the determinations of the gluon distribution. HERA jet data have 
recently been used for this purpose in the ZEUSJETS PDF fit, as described 
in Sec. 3. Clearly the combination of inclusive cross-section data and jet 
production data could lead to a simultaneous determination of as{Mz) and 
the gluon PDF. An early analysis was performed by HI ^, and this idea has 
recently been fully developed into a simultaneous fit of as{Mz) and all the 
PDFs in the ZEUSJETS-a^ fit as also described in Sec. 3. 

2. Measurements of oca from jet production data 

The method used to extract as from jet measurements is nicely illustrated 
on the left hand side of Fig. 2. A programme, like DISENT or NLOJET, 
is used to calculate NLO QCD cross-sections for jet production, for several 
fixed values of as{Mz). The parton distribution functions (PDFs) input 
to these calculations must be compatible with these as{Mz) values. (In 
practice the CTEQ4A, MRST99 and ZEUS-S variable series have been 
used). The predictions for each cross-section bin, i, are fitted to the form 

ai{as{Mz)) = Aas{Mz) + B,al{Mz), (1) 

to determine Ai and Bi so that we have a prediction for any value of 
as{Mz)- We then map the measured value of cTj and its statistical error 
onto this function to determine as{Mz) and its statistical error, as illus- 
trated in the Figure. Experimental systematic uncertainties are evaluated 
by repeating the procedure with the measurement shifted according to each 
source of systematic uncertainty and theoretical uncertainties are evaluated 
by repeating the procedure with the predictions shifted according to differ- 
ent theoretical assumptions. 

A summary of HERA as{Mz) determinations, including those presented 
at EPS05, is given on the right hand side of Fig. 2. A HERA average has 
been constructed from these results ^ , using the known experimental corre- 
lations within each experiment and assuming full correlation of theoretical 



7,2008 8:12 Proceedings Trim Size: 9in x 6in 



coopersarkar 



3 



80r 



70 



60 



50L 



NLOJET++: NLO (1+5„J 
Measured cross section o 
fit M,{MJ + BaJiMj) 
Propagation of o to as(fvlz) 



0.110 



0.115 



0.120 

«s(Mz) 



Inclusive jet cross sections in NC DIS 

ZEUS prel. (contributed paper to EPS05) 

Inclusive jet cross sections in NC DIS 

HI prel. (contributed paper to EPS05) 

Multi-jets in NC DIS 

HI prel. (contributed paper to EPS05) 

Multi-jets in NC DIS 

ZELS (DESV 05-019 - hep-ex/OS 02007) 

Jet shapes in NC DIS 

ZEUS (NucI Phys B 700 (2004) 3) 

Inclusive jet cross sections in 7p 

ZEUS (Phys Lett B 560 (2003) 7) 

Suhjet multipUcity in CC DIS 

ZEUS (Eur Phys Jour C 31 (2003) 149) 

Suhjet multiplicity in NC DIS 

ZEUS (Phys Lett B 558 (2003) 41) 

NLO QCD fit 

HI (Eur Phys J C 21 (2001) 33) 
NLO QCD fit 

ZEUS (DESY 05-050 - hep-ex/05 03274) 
NLO QCD fit 

ZEUS (Phys Rev D 67 (2003) 012007) 
Inclusive jet cross sections in NC DIS 
HI (Eur Phys J C 19 (2001) 289) 
Inclusive jet cross sections in NC DIS 
ZEUS (Phys Lett B 547 (2002) 164) 
Dijet cross sections in NC DIS 
ZEUS (Phys Lett B 507 (2001) 70) 
World average 
(S. Bethke, hep-ex/0407021) 



a{My) 



Figure 2. Left side: the method for determining as{Mz)- Right side: HERA as{Mz) 
measurements, compared to the world average 



scale uncertainties but no correlation between ZEUS and HI data. The 
result is 

as{Mz) - 0.1186 zb 0.0011(exp.) ± 0.0050(th.) 

The most accurate results come from inclusive jet production and tri- 
jet/di-jet/single jet ratios. The recent results for these processes are now 
discussed in more detail. 

ZEUS have measured inclusive jet production in DIS using S1.7pb~^ of 
data from the 98-00 running period ^. The cross-sections da/dEj^^^ where 
Ej^^ is the transverse energy of the jet in the Breit frame, are compared to 
the predictions of DISENT, for various regions, in Fig. 3. The method 
outlined above has been used to extract as (Mz). The best result is obtained 
for > 500 GeV^, when the theoretical uncertainties are minimized: 

a,{Mz) = 0.1196 ± 0.0011(stat.) ig:^!]!^ (sys.) ±^omu (th.) 

This can be compared to the previous ZEUS resuh from a similar treatment 
of data from the 96-97 running period: 

as{Mz) - 0.1212 ± 0.0017(stat.) ±"omll (sys.) i^l!'? (th.) 

Not only do the higher statistics of the newer data reduce the statistical 
error, they also allow a reduced estimate of the systematic error. The 
biggest contribution to the experimental systematic error comes from the 
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Figures. da/dExB for inclusive jet production. Leftside: ZEUS data. Right side: HI 
data 

absolute energy scale of the jets. For ZEUS this is ~ 1%, leading to a ^ 5% 
error in the jet cross-sections. The theoretical uncertainties come from 
the usual procedure of varying the renormalisation and factorisation scales 
(Pr = fip = Et) by a factor of two. Fig. 3 also illustrates similar data 
from HI from 61pb^^ of 99/00 running ^. The value of as{Mz) extracted 
from these data is 

a.iMz) = 0.1175 ± 0.0016(exp.) ±"„ °"„ll (th.) 

The method used to extract as can also be used to fit for as{< Et >) 
or as{< Q^) >■ Fig. 4 illustrates these results as a function of Q^, demon- 
strating the running of with scale, from data within a single experiment. 

Both experiments have used data from the same running periods to 
measure mutlijet production Previously the ratio of di-jet to single-jet 
cross-sections has been used to extract as- However, increased statistics 
now allows us to use tri-jet to di-jet cross-sections and these are more sen- 
sitive to tts- The advantage of the ratio technique is that many correlated 
experimental and theoretical uncertainties cancel out in ratio. The usual 
method to extract as is apphed to predictions for R^^/ R2{as{Mz))- The 
result for ZEUS is: 

as{Mz) - 0.1179 ± 0.0013(stat.) ij];™^* (sys.) ±lZf^ (th.) 
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Figure 4. Left side: illustration of running as from ZEUS analysis of inclusive jet data. 
Right side: Illustration of running as from HI analysis of tri-jet to di-jet cross-section 
ratios 

and for HI: 

a,{Mz) = 0.1175 ± 0.0017(stat.) ± 0.005(sys.) ±omlt (th.) 

As usual the major contribution to the experimental systematic error is the 
jet energy scale uncertainty, and for the theoretical uncertainty it is the 
remormalisation and factorisation scale uncertainty. Fig. 4 also illustrates 
the running of ag from the HI analysis of multijet production. 

3. Use of jet production data to determine the gluon and 
otsiMz) in PDF fits 

Both HERA collaborations have recently made PDF fits using exclusively 
their own data The advantages of a fit within one experiment are that 
the systematic uncertainties are well understood, and the advantages of 
a HERA only fit are that these data have the most extensive kinematic 
coverage of any DIS data, on a pure proton target, hence without the un- 
certainties of heavy target corrections, and without the need to make strong 
isospin assumptions. However, such HERA only analyses have only been 
possible since the 98/00 running extended HERA kinematic coverage to 
high , where the contributions of and Z exchange become impor- 
tant. This allows us to gain information on the valence d and u quarks at 
high X, from the e'^p and e~p charged current cross-sections, respectively. 
Furthermore, the difference between the e^p and e~p high-Q^ neutral cur- 
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rent cross-sections allows us to extract the valence structure function xF^ 
for all X. Valence information was the missing piece in a HERA only anal- 
ysis, since HERA data have provided the major contributions to the mea- 
surement of the low-a; sea and gluon distributions, through the structure 
function, F2, and its scaling violations, dF2 / dln{Q'^) , respectively. It has 
been a limitation of all PDF fits, including the MRST and CTEQ global fits 
to world data, that information on the high-x gluon is lacking. The global 
fits have used Tevatron jet data to remedy this. However, the HERA jet 
data have the advantage of much smaller experimental uncertainties and 
do not suffer from possible uncertainties due to new physics. 

The ZEUSJETS PDF fit uses all the inclusive double differential cross- 
section data from HERA-I running (94-00) and adds the DIS inclusive jet 
production data from ^ and the di-jet data from photoproduction ^° from 
96-97 running, to gain information on the gluon. As remarked earlier, the 
jet cross-sections depend on the gluon and quark PDFs and on a,, and 
this can be used to break the strong coupling between and the gluon 
PDF which plagues the indirect extractions from the scaling violations of 
inclusive data. The fit takes full account of the correlated systematic un- 
certainties within and between the data sets and these are propagated into 
the PDF uncertainties using the OFFSET method. 

A remark is in order on the use of photoproduced di-jet data. It is well 
known that at leading order photoproduction can proceed via direct and 
resolved processes. In the former, the photon interacts directly as a point- 
like particle, whereas, in the latter, the photon has decomposed into quarks 
and it is one of these quarks which initiates the interaction. The resolved 
cross-sections will have some sensitivity to the photon PDF. We avoid this 
by restricting the analysis to direct enriched cross-sections by the cut 

« = i:iEl'*'exp{-rf'**)/2yEe > 0.75 

Of course at NLO these distinctions cannot be made precisely, but this cut 
still serves to minimise the influence of the choice of photon PDF. 

A new technique had to be developed to include NLO jet cross-sections 
in a PDF fit, since the computation of such cross-sections is very CPU 
intensive. Hence the NLO programmes were used to compute LO and NLO 
weights, fT, which are independent of as and the PDFs, and are obtained by 
integrating the corresponding partonic hard cross sections in bins of ^ (the 
proton momentum fraction carried by the incoming parton), fip and fiji. 
The NLO QCD cross sections, for each measured bin, were then obtained 
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Figure 5. PDFs from the ZEUSJETS PDF fit. Left side: u and d valenee PDFs. Right 
side: sea and gluon PDFs 

by folding these weights with the PDFs and according to the formula 

n a ij,k 

where the three sums run over the order n in as , the flavour a of the incom- 
ing parton, and the indices {i,j,k) of the ^, /^f and fin bins, respectively. 
The PDF, fa, and as were evaluated at the mean values (^), (fip) and 
(fJ-B,) of the variables ^, fj-p and /ifl in each {i,j,k) bin. This procedure 
reproduces the NLO predictions to better than 0.5%. 

The ZEUSJETS PDF fit gives a very good to aU of the input data 
sets, simultaneously describing inclusive cross-sections and jet production 
data, thus providing a compelling demonstration of QCD factorisation. The 
extracted PDFs are shown for various pig. 5. The ZEUSJETS PDFs 

are in good agreement with world PDF extractions and are shown compared 
to the HI PDFs in Fig. 6. This right hand side of this figure illustrates how 
the input of the jet data has improved the uncertainty on the gluon PDF, 
by about a factor of two, in the region 0.1 < a; < 0.4. 

Most world PDF fits are performed assuming a fixed value of as{Mz) 
equal to the world average (asiMz) = 0.118 for the ZEUSJETS fit de- 
scribed above). This is because the coupling between the gluon density 
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Figure 6. Left side: ZEUSJETS PDFs compared to HI PDFs 2000. Right side: the 
fractional uncertainties on the gluon PDF are shown before and after the jet data are 
included in the ZEUSJETS PDF fit. 
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Figure 7. Left side: profile for as{Mz) with and without jet data. Right side: the 
ZEUS JETS— as PDFs with the uncertainty from variation of aa{Mz) included. 



and the value of as{Mz) in the DGLAP equations does not allow a very 
accurate measurement of as{Mz)- However, with the added information 
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Figure 8. Projected PDF improvements after HERA-II running. 



on the gluon which comes from the jet data, we can make an accurate ex- 
traction by allowing as{Mz) to be a free parameter simultaneously with 
all the PDF parameters, in the ZEUSJETS-a^ fit: 

as{Mz) = 0.1183±0.0007(stat.)±0.0027(sys.)±0.0008(model)±0.005(th.) 

where the systematic error contains all the correlated experimental system- 
atics of the contributing data sets. The model error includes uncertainties 
due to changes in the choice of the proton PDF parametrisation at the 
starting scale, and in the choice of the photon PDF. The theoretical error 
comes from variation of the choice of the renormalisation and factorisation 
scales as usual. Fig. 7 illustrates the contribution of the jet data to the ac- 
curacy of the as extraction by comparing the profile, with and without 
the jet data. This is the first extraction from HERA data alone because 
the extractions described in Sec. 2 all require input of PDFs which were 
extracted using other data. In the ZEUSJETS fit we are able to vary the 
as(Mz) used for predicting the jet cross-sections, simultaneously with vary- 
ing the PDFs. A consequence of this is that the uncertainties on the PDFs 
due to the uncertainty on as{Mz) is much reduced in the ZEUSJETS-as 
fit, as also illustrated in Fig. 7. 
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4. Outlook 

We can look forward to improving these PDF fits with more jet data in the 
near future, since the jet data included in the present fit came only from 
the 96-97 running. Jet data from the 98-00 running are already becoming 
available, as described in Sec. 2, where they were used for preliminary 
extractions Combination of ZEUS and HI data will also yield greater 
accuracy. In the mid-term future we look forward to improving these mea- 
surements with HERA-II data. A study of possible improvements has 
assumed the achievable goals of a luminosity of 350p6~^ for both e"*" and e~ 
running, and BOOpb^ of jet data. The extra statistics of high data will 
yield improved high- a: valence information, and an extraction of xF^ will 
extend the valence information down to a; ~ 0.01. The extra statistics will 
also improve the high-.T sea extraction, whereas the jet data will further 
improve the high-x gluon. Fig. 8 shows the expected level of improvement 
in PDF uncertainty. 
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